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We analyze behavior of correlated electrons described by Hubbard-like models at 
intermediate and strong coupling. We show that with increasing interaction a pole in a 
generic two-particle Green function is approached. The pole signals metal-insulator tran- 
sition at half filling and gives rise to a new vanishing "Kondo" scale causing breakdown 
of weak-coupling perturbation theory. To describe the critical behavior at the metal- 
insulator transition a novel, self-consistent diagrammatic technique with two-particle 
Green functions is developed. The theory is based on the linked-cluster expansion for 
the thermodynamic potential with electron-electron interaction as propagator. Parquet 
diagrams with a generating functional are derived. Numerical instabilities due to the 
metal-insulator transition are demonstrated on simplifications of the parquet algebra with 
ring and ladder series only. A stable numerical solution in the critical region is reached by 
factorization of singular terms via a low-frequency expansion in the vertex function. We 
stress the necessity for dynamical vertex renormalizations, missing in the simple approxi- 
mations, in order to describe the critical, strong-coupling behavior correctly. We propose 
a simplification of the full parquet approximation by keeping only most divergent terms 
in the asymptotic strong-coupling region. A qualitatively new, feasible approximation 
suitable for the description of a transition from weak to strong coupling is obtained. 



I. INTRODUCTION 

Tight-binding models of correlated electrons are expected to provide description of thermodynamic and 
spectral properties of the underlying system in weak as well as in strong coupling regimes. We have at our 
disposal relatively reliable techniques to describe the two extreme limits of weak and strong couplings in 
the archetypal Anderson impurity and Hubbard lattice models. The weak-coupling regime is governed by a 
Hartree-Fock mean field with dynamical fluctuations captured by Fermi-liquid theory. Extended systems 
at low temperatures are Pauli paramagnets with smeared out local magnetic moments. For bipartite 
lattices antiferromagnctic long-range order sets on at half filling at arbitrarily small interaction. In the 
strong-coupling regime the Hubbard model at half filling maps onto a Heisenberg antiferromagnet with 
pronounced local magnetic moments and the Curie- Weiss law for the staggered magnetic susceptibility, at 
least at the mean-field level. The spectral structure is dominated by separated lower and upper Hubbard 
bands and the strongly correlated system seems to be insulating even in the paramagnetic phase. 

In recent years importance of the regime, where the effective Coulomb repulsion is comparable with 
the kinetic energy and hence neither very weak nor very strong, has significantly increased. Unfortunately 
there is up to now no adequate method for the description of this transition regime between the weak 
and strong-coupling limits. At intermediate coupling dynamical fluctuations control the low-temperature 
physics of interacting electrons and neither weak-coupling nor atomic-like perturbation theories are ade- 
quate. In this nonperturbative regime we expect the Kondo regime to set on in the Anderson impurity 
model and eventually breakdown of the Fermi-liquid regime and metal-insulator transition in the 
Hubbard model §. 

New progress in the construction of dynamical mean-fleld theories via high spatial lattice dimensions 
pushed forward our possibilities to investigate the transition region between weak and strong coupling 
1^. Although the spatial fluctuations are reduced in high dimensions, there is no restriction on quan- 
tum dynamical fluctuations important at intermediate coupling. Particularly the quantum Monte-Carlo 
technique in d — 00 revealed a number of new features of the (disordered) Hubbard model at finite tem- 
peratures 0-0. However, the quantum Monte-Carlo is restricted only to relatively high temperatures. 
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That is why an analytic method based on second-order perturbation theory, caUed IPX, was used to study 
the transition from weak to strong couphng at zero temperature ||^. IPX at half filling reproduces well the 
finite-temperature Monte-Carlo data and contains the weak-coupling (up to U^) and atomic solutions as 
exact limits. Contrary to Hubbard-III-like theories it reproduces Fermi liquid at weak coupling. A metal- 
insulator transition was found at Uc ~ G/Tru, where u is the DOS at the Fermi level ||9|-[ll|]. Although IPX 
is an analytic theory, its solution can be reached only via iterations that depend on the initial value of 
the self-energy. When we start with £'^"^(2;) = the iterations converge for weak interaction to a metallic 
solution. If we start with E(o)(z) = Eat(^) = U^/4:Z we end up with an insulating solution for sufficiently 
strong interaction. Xhe IPX scheme fails to converge close to the metal-insulator critical point. Another 
reasoning had to be used to demonstrate the existence and character of the metal-insulator transition 

At present IPX seems to be the only analytical-numerical scheme interpolating between Fermi liquid at 
weak coupling and paramagnetic insulator at strong coupling in a frustrated system, where the antiferro- 
magnetic order is suppressed. However, the solution is numerically unstable in the most interesting region 
around the critical interaction. Xhe quantitative details of the critical behavior at the metal-insulator 
transition and the way Fermi-liquid theory breaks down remain unaccessible. IPX is not a conserving 
theory and the two-particle functions, important for the description of inter-particle correlations, are 
not reliable in the IPX scheme. A systematic, consistent, and controllable theory delivering quantitative 
description of the critical behavior at the metal-insulator transition is still missing. 

Xhe aim of this paper is to fill up this gap and to develop a systematic theory interpolating between 
Fermi-liquid and local-moment regimes suitable for the description of the critical behavior at the metal- 
insulator transition of the Hubbard model or the Kondo behavior in the single-impurity Anderson model. 
We start with weak-coupling, self-consistent perturbation theory and propose a scheme how to extrapolate 
it to strong coupling. Contrary to IPX we concentrate on the two-particle Green functions showing 
divergence at the metal-insulator transition. We develop a systematic renormalized perturbation theory 
based on parquet-type diagrams with nontrivial renormalizations of two-particle functions. Especially 
we demonstrate the necessity for dynamical vertex renormalization whenever we approach a critical 
point. A theory with sufficient vertex renormalization is derived to describe the critical region of the 
metal-semimetal transition with condensation of bound electron-hole pairs. It is a simplified version of 
the parquet algebra keeping only most divergent contributions to the two-particle scattering function 
diverging at the critical point. It strongly renormalizes scatterings within singlet electron- hole pairs 
and leaves inter-pair and triplet correlations unrenormalized. Xo stabilize the numerical solution in the 
critical region the two-particle pole is factorized in the spirit of the renormalization group. Xhe singular 
contributions can then be integrated analytically. Xhe resultant new equations allow for a numerically 
stable solution even in the critical region of the condensation of bound electron- hole pairs. Such a method 
of numerical evaluation is asymptotically exact in the critical region of effective impurity models within 
the chosen approximation and makes tractable otherwise inaccessible diagram sums renormalizing the 
coupling constant far beyond simple ladder or ring series. 

Xhe layout of the paper is as follows. Section II analyzes the reasons why numerical solutions of 
sample models break down when the electron-electron interaction becomes strong. A systematic pertur- 
bation theory for intermediate and strong coupling using two-particle functions is expounded in Sec. HI. 
Simple approximations of ring and ladder diagrams are rederived as examples in Sec. IV. Xhese simple 
approximations are used in Sec. V to demonstrate the existence of poles in two-particle functions im- 
peding direct numerical continuation of weak-coupling solutions to the strong-coupling limit. Section VI 
demonstrates the necessity for vertex renormalizations in self-consistent theories around critical points in 
order to comply with a Ward identity binding charge renormalizations (vertices) to mass renormalizations 
(self-energy). A feasible approximation within parquet algebra with self-consistent vertex renormaliza- 
tions suitable for the description of the transition from weak to strong coupling is derived in Sec. VII. 
Last Section VIII brings conclusions. 

II. CORRELATED ELECTRONS AT INTERMEDIATE COUPLING 

Xo study correlated electron systems at intermediate coupling quantitatively we will basically work 
parallelly with two sample models: the single-impurity Anderson and the lattice Hubbard models. Xhe 
former one will serve as a test case for comparison of developed approximations with the exact Bethe- 
ansatz or numerical renormalization-group solution. Xhe latter model will explicitly be studied only at 
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the mean- field level, i.e. in infinite spatial dimensions where spatial fluctuations are suppressed |l3|| . This 
enables us to cover both the model situations within one theoretical framework. 
We start with the underlying Hamiltonians 

HsiAM = E ^(k)cLck. + ef J2 fU^ + Ul\f^ fU, + E [^kcL/. + /^^J (la) 

k,cr (J k.cr 

for the single-impurity Anderson model and with 

k,a" i 

for the Hubbard model. In the former case we integrate out the conduction electrons and transform the 
model onto a single-site theory with dynamical localized electrons. The interaction dependent partition 
sum has form of a local functional integral 

Zloc = j VfVf* cxp 1^ f;nQ^,\Un - U J^^ drf^ {T)f^ (t)/*^ (t)/^ (t) | (2a) 



where 



g-,l = iLOn + ^Ia-ef- V^G4iLUn) , Gi'>\z) = J dep{e)[z + - e]-^ (2b) 

is the effective dynamical propagator of the localized electrons and fj,^ = ii + ah . The variables /, /* are 
Grassmann ones and p is the density of states of the conduction electrons. 

The Hubbard model in infinite dimensions, i.e. at the mean-field level can be mapped onto an 
impurity problem self-consistently embedded in an averaged medium jlj,^ . The partition sum can then 
be transformed onto a local functional integral identical with ( |la|) where only the effective dynamical 
propagator is defined as 

g^jl = G^,] + Eo-n , Gan = J dep{e)[iuJn + fla - ^an - f]"^ (2c) 

with So-n as the self-energy to be determined from an equation 



G.„ = y^ / I?/2?r/:„/anexp|^/;„C;-„V.„-{/ f dTf*{T)f^{T)fl{T)f^{T) 
'-hoc J I „ Jo 



(2d) 



The physics of both the sample systems is determined completely from the effective local but dynamical 
partition sum 7-hoc- Here the unperturbed propagator is Q and the only dynamical variable is the 
imaginary time r e (0,/3) or Matsubara frequency za;„. 

Although the mathematical representation for the impurity Anderson and the mean-field Hubbard 
partition sum is formally identical, physical behavior of these models at intermediate coupling is expected 
to be different. The impurity model remains local Fermi liquid up to infinite interaction strength [ p^ . 
At half filling a quasiparticle peak forms around the Fermi energy (Ep = 0) in the density of states of 
the local dynamical electrons the width of which is determined by a Kondo (dimensionless) scale 

A = exp{-7rC/i^} (3) 

and a typical energy of the model (e.g. the bandwidth of the conduction electrons). It is the existence 
of this Kondo scale incomparably smaller than other scales of the underlying Hamiltonian that makes 
perturbation expansion fail at intermediate coupling. Except for asymptotically exact solutions in the 
low-energy limit such as Bethe-ansatz or numerical renormalization group [^ , there is no other 
reliable approximate scheme being able to reproduce the Kondo strong-coupling asymptotics (^) . 

Recent studies on the Hubbard model in d = cxd show that the Kondo scale of the impurity model 
should vanish at a critical interaction (JJc ~ 67r/y) due to the additive self-consistency between the "un- 
perturbed" propagator Q and the self-energy p-p!^. The electronic system undergoes a metal- insulator 
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transition. Fermi-liquid theory and weak-coupling perturbation expansions break down at this transi- 
tion. Breakdown of perturbation theory at the Mott metal-insulator transition is not caused by the 
exponentially small Kondo scale but by the divergence of the effective mass determining criticality of 
the metal-insulator transition at zero temperature. Effective mass is connected with the real part of 
the self-energy m* /m cx 1 — E' where E' = 9S(a;)/9a;|i^=o- Divergence of the effective mass also causes 
vanishing of an energy scale. This energy scale is determined by frequencies uj± at which the real part of 
the self-energy reaches its extremum closest to the Fermi energy. These frequencies determine an energy 
interval around the Fermi level within which Fermi liquid and the quasiparticle picture is applicable. 
Fig. 1 shows ReS(cL') calculated within self-consistent ring diagrams. It is clearly seen that uj± with 
increasing interaction strength. Vanishing of the difference between oj+ > and uj- < and arising 
of a jump in Rel](a;) then lead to numerical instabilities and breakdown of iteration schemes to self- 
consistent approximations at intermediate coupling, disregarding whether they can or cannot describe 
the metal- insulator transition |19|. 

At intermediate coupling the two investigated models have a vanishing small scale in common. It 
causes numerical instabilities in self-consistent approximations and breakdown of iteration schemes [ pO| . 
It hinders analytic continuation of weak-coupling approximations to the strong-coupling limit. To suc- 
ceed in quantitative description of intermediate coupling it is necessary to treat the vanishing "Kondo" 
scale separately from the other scales determined by the underlying Hamiltonians. To do that we must 
understand the origin of this scale. We will demonstrate in this paper that the vanishing scale reflects 
criticality of a two-particle (actually electron-hole) correlation function. It means that a two-particle 
Green function at intermediate coupling develops a pole at the Fermi energy. It is known from the early 
papers on the Kondo problem that the Kondo scale or the Kondo temperature is connected with an RPA 
pole in the magnetic susceptibility It was pointed out recently Q that also at the metal-insulator 
transition the local, dynamical susceptibility shows divergence at the Fermi energy. In order to get the 
vanishing scale at intermediate coupling under control a systematic renormalized perturbation theory 
based on approximations for two-particle Green functions is to be used. 



III. SYSTEMATIC SELF-CONSISTENT APPROXIMATIONS FOR TWO-PARTICLE GREEN 

FUNCTIONS 

To develop systematic approximations for one-particle Green functions we use Dyson's equation en- 
abling us to consider explicitly only one-particle irreducible diagrams contributing to the self-energy. The 
situation is more complicated at the level of two-particle functions. Although we can use a Bethe-Salpeter 
equation to extract two-particle irreducible diagrams, we have three possibilities to do it. We have three 
channels for multiple two-particle scatterings defining three types of two-particle irreducibility. They 
are the electron- hole (e-h), electron-electron (e-e), and interaction (U) channels schematically drawn in 
Fig. 2. Different channels mean different rearrangements of perturbation expansion. The chosen irre- 
ducible diagrams are summed first and the remaining reducible ones are summed via a Bethe-Salpeter 
equation in the end. If perturbation series converges all rearrangements must lead to the same result. 

When studying two-particle functions we need to know the one-particle self-energy. Hence perturba- 
tion theory for two-particle functions cannot be developed without a parallel expansion for the self-energy. 
To approximate the one and two-particle functions simultaneously we use perturbation expansion with 
two-particle functions as a means for developing controllable comprehensive approximations for a gener- 
ating thermodynamic potential. All the physical quantitities will then be determined consistently from 
the generating functional via functional derivatives. 

As a first step we decide which is the relevant two-particle function for our purpose to describe 
intermediate and strong coupling. It should be a function approaching a pole with increasing interaction. 
Such a function must be related to the interaction (two-particle) part of the underlying Hamiltonian. 
The natural generalization of the Hubbard interaction is the following two-particle function 

Cij(r) = (nii(r)Hjx(O)) ~ {m^{r)){ni^m (4) 

measuring correlations between spin up and spin down densities at different times and different lattice 
sites. The brackets denote thermal average. Fourier transform C{q,ii'm) of the above function is a 
two-particle function approaching a pole when the Hubbard interaction increases. Note that Vm are 
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bosonic Matsubara frequencies. Due to the space and time translational invariance of the Hamiltonian, 
ImC(q, 0) ~ and the real part has definite sign, namely 



C(q,0) (X -UXM)X,{ci) < with XAci) - AA-i Ek ^^^^'^^Ikjlfffq')'^^'"'^ ■ (5) 

This relation is exactly fulfilled at weak-coupling and can qualitatively be broken (change of the sign) 
only if the function C(q, goes thru a pole. The pole would indicate breakdown of Fermi liquid. In 
the weak-coupling, Fermi-liquid regime we can write 

C(q,0) = ^[K(q,0)-x(q,0)] <0 (6) 

where ^(q, ii^m) is the dynamical compressibility and x(Q7*^m) the dynamical magnetic susceptibility. 
Eq. (^) states that magnetic fluctuations are stronger than charge fluctuations in the models with the 
Hubbard interaction. Divergence in the function C(q, 0) induces divergence in the magnetic susceptibility 
and hence is indication of a magnetic long-range order. 

Function C(q, ivm) will be approximated by a selection of classes of Feynman diagrams. Consistency 
and conservation laws demand that a two-particle approximation be in concord with an approximation 
on the one-particle self-energy. According to Baym and Kadanoff |^,^ a conserving theory must work 
with fully renormalized one-particle propagators. Hence only fully self-consistent approximations are 
acceptable at the level of two-particle functions. To construct consistent approximation schemes we have 
to build up the generating functional ^[U, G] using the function C(q, Wm)- To do that we apply a linked- 
cluster expansion recently proposed by the author and J. Schlipf |20[| . Its salient feature is that the the 
interaction line (non-relativistic photons) is used for edges (diagram bonds or propagators) and loops 
of electron Green functions serve as unperturbed vertices. An analogy between the classical free-energy 
functional of the Ising model VF[J, /i] and the quantum grand potential ri[t/, G^"^] is thereby made, i.e. 
J < — > U and h < — > G^"-*. We use this analogy to construct the generating functional $[C/, G] from the 
two-particle Green function C(q, Wm). 

We represent the grand potential of the Hubbard model in the homogeneous phase in the following 
form 

S, G] = -Un^ni - ^ g»",iO+ ^^j^ [i^j^ + fia - e(k) - Un^„ - Ecr(k, iujn)] 

cjn,k 

+G,(k, iw„)S,(k, iuj^)} + $[[/, G] , (7) 

where n^, 5]o-(k, iw„), Go-(k, iujn) are variational variables (functions). We now introduce a small pertur- 
bation U U+6U (q, ivm)- The function C(q, Wm) is defined as a variational derivative of the generating 
functional 



(5C/(q, Wm) 

Linked-cluster theorem can be used for the inverse transformation 



(8) 

su=o 



qm ^ 



(9) 



Two-particle function C(/7A|q, ivm) is understood as a functional of the full renormalized electron prop- 
agator G(k, iuJn) with the self-energy determined at the interaction strength U . 

Function C{U\(\,iUm) does not directly obey an equation of motion. It is connected with a general 
two-particle Green function /C-f|(C/|fc', A:"; g) determined from a set of coupled Bethe-Salpcter equations. 
To simplify lengthy expressions we sometimes use a four-vector notation k = (k, iuJn) and q = (q, ivm) 
for fermionic and bosonic variables, respectively. We can write 

C(C/|q,zi/™) = -i^ ^ G^{k')G^{k' + q)G^{k")G^(k" + q)]C^^{U\k\k"-q). (10) 

k' ^k" 

Functions Kaa' (k' , k"; q) can be represented as a sum of contributions from different two-particle channels 
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K,,. (k, k'; q) = (fc, k'-q)+ IC^,, {k, k'-q) + ICf,, (k, k'; q) + /C^ - {k, k'; q) (11) 

where I^a' contains diagrams irreducible in all the channels while K,"^, is a sum of diagrams reducible in 
the particular channel a. Each of the reducible two-particle functions obeys a Bethe-Salpeter equation 
mixing all the channels and eventually spins. The equations in the "horizontal" electron-hole and electron- 
electron channels read 

Kt, {k, k';q)^-^Y. + ik' + q') 

q' 

lC>{k + q',k' + q':q-q'), (12a) 
/C- , (fc, fc'; g) ^ - -1^ ^ A- - (fc, A; + fc' + g - k" ■ k" - k)G, {k")G,, {k + k' + q - k") 

k" 

lCaa'{k",k'-k + q-k"). (12b) 
The reducible functions in the "vertical" or the interaction channel are determined from equations 

K'i,,{k,k'-q)^-^ J2 A';^^4k,k";q)G,.{k")G,„{k" + q)IC„„,,{k",k';q). (12c) 

k" a" 

Here A^^, (fe, k'; q) are sums of irreducible diagrams in the a-th channel. Equations ( [T^ ) have the structure 
of parquet diagrams ||2j,|2^ . To complete the parquet algebra and to close the equations we must add a 
definition for the irreducible functions from each channel. When no further approximations are used the 
parquet equations are completed with a relation between reducible, JC", and irreducible, A" functions 

A^,, {k\ k"- q') = K„„, (fc', fc"; q) - K^,, (fc', fc"; q) . (13) 

The generating functional <&[[/, G] is now fully determined from the completely irreducible two-particle 
vertex functions /g-o-' (C/A|fc', A:"; g), Bethe-Salpeter equations ([III)- ([l3|), and the linked-cluster theorem 
(|^), with (p^). The self-energy is determined from the saddle-point of the generating functional VL with 
respect to the fermion propagator G. Note that the self-energy here measures only dynamical fluctuations 
beyond the Hartree approximation. Even the simplest approximations on the scattering matrix /Cf ^ lead 
to nontrivial complex structure of two-particle functions causing the existence of poles at strong coupling. 



IV. SIMPLE APPROXIMATIONS 



As examples of the general construction we derive a few simple approximations where we can get 
rid of the integral over the interaction strength and calculate all the quantities at the same value of the 
coupling constant U . In all these approximations the completely irreducible vertex /f| is assumed to be 
the bare interaction U . 



A. Ring diagrams 

We neglect contributions to the two-particle scattering matrix /Cf ^ coming from the electron-electron 
and electron-hole channels and consider explicitly only the interaction channel. It means that no multiple 
scatterings are allowed. Further on we neglect the irreducible functions for scatterings of quasiparticles 
with the same spin. The two-particle function /Cf| in this approximation then reads 

The integral over the interaction strength can simply be performed leading to a generating functional 

^^"^'[C/, G](q, ivm) = ^ E e"""^ 1^ [1 - U^^M^ ^Um)X^{cl, lu^)] (15a) 

qm 
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where 



X„(q, w„) = ■^^Ga{k,iujn)Ga{k + C[,i{uJn + lym)) (15b) 

kn 

is a bubble of two full one-electron propagators. The self-energy can easily be derived in form of an 
integral equation 



2pAf 



qm 



We have a complete set of integral equations the solution of which fully determines thermodynamic as well 
as spectral properties of the self-consistent ring-diagram approximation. If we reduce this approximation 
to infinite dimensions, the full extended propagator is replaced by its local element, we reveal a mean-field 
version of the ring diagrams recently studied numerically At the mean-field level the numerics 

breaks down at intermediate filling {U « G/tti/) when the pole in C{z) is approached We expect the 
same for the finite-dimensional version, but the pole will appear at a specific wave vector qg. 

B. Ladder diagrams 

Next step in developing approximations onto two-particle scattering matrix /Cf ^ is to consider ladders 
of either multiple electron-hole or electron-electron scatterings by neglecting the contributions from the 
interaction channel. In these ladder approximations we neglect screening of the interaction due to polar- 
ization bubbles of electron-hole pairs. Summing the geometric series of the Bethe-Salpeter equations we 
obtain for the electron-hole channel 

/Cf|^^(f/|fc, k'- q) = ^— ^ ^ — (16a) 

l + ^Efe"GT(fc-fc' + A:")Gx(fc") 

and for the electron-electron channel 

Klt'^{U\k, k'; q) = r-j ^ ^ ^ • (16b) 

^ ' ' '"^^ l+v_^Y.k"G^{k + k'-k")G^{k") ^ ^ 

We can again perform the integration in the linked-cluster theorem (^ and end up with 

cE>^^^[C/, G] = ^e*-"-o+ {UX{ci,iVm) - ln[l + C/X(q, ii/^)]} (17a) 

qm 

for the electron-hole channel and analogously with 

qm 



for the electron-electron channel. Here we denoted 

X{(i,Wm) = ^^G|(q + k,itj„+„)Gx(k,iw„), (17c) 



k7i 



y (q, = ^ X! (q - iuJni-n)Gl (k, iujn) . (17d) 



kn 



Note that generates random-phase approximation with rcnormalized propagators, i.e. its renor- 

malized version as introduced by Suhl in the study of the single- impurity Anderson model pTf . The gen- 
erating functional with multiple electron-electron scatterings describes a rcnormalized T-matrix 
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approximation studied by Baym and Kadanoff ||2^. The self-energies, respective their corrections to the 
Hartree approximation, for ladder diagrams read 

Er^(k, ^c.„) = ^ G_.(k - q, ^u.^_„:)^S^p}L- , (17e) 

/37V ^ 1 + C/X(q,2z/™) 

f-' 1 + C/y(q,zj/„,.) 



f3M 



The ring and ladder diagrams constitute simplest approximations at the level of two-particle functions 
enabling for explicit analytic form of the generating functional These approximations are fully self- 
consistent and contain geometric series of two-particle bubbles. The ring diagrams use bubbles of electron- 
hole pairs of the same spin, while the RPA and TMA involve bubbles of quasiparticles with opposite spins. 
Due to the full self-consistency in these simple approximations it is in principle possible to find a solution 
in the whole range of the interaction strength. It is, however, necessary to stress that these approximations 
can be expected reliable only at weak coupling, since only diagrams to second order are treated exactly 
here. We use these approximations as samples to demonstrate breakdown of numerical iteration schemes 
at intermediate coupling and the necessity for a separate treatment of the new small scale connected with 
zero of the denominator of the two-particle functions from the self-energy equations. 



V. TWO-PARTICLE CRITICALITY 



To simplify the analysis we leave the general formulation and resort to a mean-field description. We 
suppress the momentum convolutions and come from the full momentum-dependent Green functions to 
their local parts. Further on we choose the spin-symmetric solution at zero temperature for the half-filled 
band. These assumptions lead to significant simplifications of the resultant equations and enable to trace 
down breakdown of perturbation theory at strong coupling of the lattice models. 

We first analytically continue the equation for the self-energy to real frequencies and take explicit 
advantage of the electron-hole symmetry. Then the generic equation for the self-energy can be written 
as @ 

= ^ / {lmG(c^;) [C{lo' - z) - C{lu' + z)] + ImC(L^V) [G{lu' + z) ~ G{u;' - z)]} (18a) 



where 



The two-particle bubble X has a representation 



XiO = - [ — ImG(c^V) [Gicj' + 0+ Giu;' - C)] 



(18c) 



Here z and stand for complex energies in fermionic and bosonic functions, respectively. 

Troubles with the numerical solution of these equations begin when the two-particle function C(C) 
approaches a pole at C = 0. In analogy to the stability condition (^) this function must fulfill 

G(0) < (19) 

in order to keep the effective mass of the interacting electrons positive. When inequality (^9|) is broken, 
even by intermediate iterations, the numerical solution fails to converge or a solution on an unphysical 
sheet is reached. It is easy to verify that A"(0) < and hence only the RPA and ring diagrams run into 
troubles at intermediate coupling. The TMA is harmless and can be extrapolated to strong coupling 
without meeting numerical problems ||28|| . One can expect this, since due to the Coulomb repulsion the 
electron-hole pairs tend to bind. Multiple scatterings of electrons on electrons do not contribute to this 
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process. Only the ring and RPA diagrams are relevant for creation of electron-hole bound states and 
divergence of the effective electron mass at the metal-insulator transition. 

The dynamical two-particle function C, the local element of the full correlation function C, defines at 
intermediate coupling a new dimensionless scale 

A = f + UX{0) (20) 

when its pole is approached. To determine the actual strong-coupling asymptotics of the RPA and ring 
diagrams it is necessary to calculate the behavior of the small scale A as a function of the interaction 
strength as accurate as possible. On the one hand the existence of this small scale impedes continuation 
of weak-coupling solutions to strong coupling. On the other hand vanishing of this scale at the critical 
point of the metal-insulator transition enables us to pick up divergent contributions in the critical region 
analytically. Only the low-frequency behavior of the two-particle function C{lu) is relevant in leading 
order of diverging effective mass. We can replace the full frequency dependence of the denominator of the 
function C(w) by its low- frequency asymptotics. Such an approach to derivation of the strong-coupling 
asymptotics was used by Hamann pO| ] in the case of Suhl's renormalized RPA and recently by the author 
191 for the ring diagrams. This treatment of the the singular function C{uj) is asymptotically exact at 
the critical point with nonintegrable singularity. 

The existing asymptotic expansions for real frequencies can be generalized to the complex plane in 
such a way that the analytic structure of the resultant self-energy is not violated. No spurious poles are 
generated in the lower and upper half-planes of complex energies. This property is essential not to break 
the sum rules. We hence make the following low-frequency ansatz for the denominator of C(C) vanishing 
at the critical point 

1 + UX{C) = A - isign(ImC)C/A:'C . (21) 

We introduced a new parameter X' as the only part of the function X{(^) relevant for the critical behavior 
of the solution. It is easy to see that for the simple approximations 

X' = Tiv'^ >Q (22) 

where v = p(0) remains unrenormalized in the Fermi-liquid regime. Hence the two-particle vertex does 
not dynamically renormalize in the low- frequency limit for simple ladder and ring diagrams. 

We use ansatz ( ^l| ) and simplify the equation for the self-energy by keeping only leading-order terms 
due to the singularity of C(C), i.e. in | In A|. As a result we obtain two regions for the self-energy. For 
finite frequencies, \z\ ^ /S./UX' the solution reads 

= is^^^^(z) = -l^lnA. (23a) 
In the limit of small frequencies, \z\ « /S./UX' the self-energy has another representation 

E«».(.) = is«'^^(z) = -gjl In A + ln[A - ^r^UX'z] . (23b) 

Formulas ( ^3| ) hold for any external (complex) frequency z = uj + irj with = signr;. The solution has 
correct Fermi-liquid asymptotics for small as well as for large frequencies. Moreover, the Herglotz analytic 
properties of the self-energy are preserved. 

To complete the equation for the self-energy we need to determine the parameter A. We use a 
representation X{Q) = 2 j^^duj' p{ui')^eG{uj') that can easily be evaluated from the derived asymptotic 
form of the self-energy. Only frequencies > oj > —\/\ In Aj/aTrA"' are relevant, where a = 1, 2 for RPA 



and ring diagrams, respectively. Using (23a) we finally obtain 



X«"'9(0) = V2X^^^(0) = J . (24) 



Note that the low- frequency asymptotics of the self-energy, (23b), is irrelevant for evaluation of the 



integral A"(0). It, however, determi nes t he Fermi-liquid, quasiparticle characteristics of the solution. E.g. 



the effective mass comes out from (|23b[ ) as 
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We know from earlier studies that Fermi liquid gets less and less thermodynamically important with 
increasing interaction, since the width of the energy interval around the Fermi level, within which Fermi- 
liquid (quasiparticle) picture holds shrinks to zero as A|lnA|. Non-Fermi-liquid contributions, i.e. the 
frequencies beyond the first extremum in the real part of the self energy overtake the control over macro- 



scopic, integral quantities [[19 31 1. 

From the definition of the parameter A, (po|), and from (|2^ ) we obtain an asymptotic form for this 
two-particle scale 

A = exp{-a7r2i/2;72} (26) 

the result derived for the first time by Hamann for the single-impurity Anderson model ]30| ] . Our deriva- 
tion of the strong-coupling asymptotics holds for the Hubbard model in d = oo as well. It means that 
there is no difference in leading singular order between impurity and mean-field lattice models at the 
two-particle critical point. 

The exponential asymptotics ( p6| ) lies one order false compared with the exact Kondo asymptotics 
(|^). The failure to reproduce the correct Kondo scale disqualifies the simple ring and ladder diagrams 
as appropriate for intermediate and strong couplings. It is evident from the parallel comparison of two 
topologically inequivalent approximations (rings and ladders) leading to the same low-frequency divergence 
that in each approximation important diagrams contributing the same divergence were neglected. We 
hence have to go beyond the simple ladder and ring diagrams in the search for solutions aiming at a 
reliable description of the transition from weak to strong coupling. 



VI. NECESSITY FOR RENORMALIZATION OF THE INTERACTION STRENGTH 



Already early studies following the pioneering papers of Suhl and co-workers and Hamann, [ p4|j32| | 
indicated the necessity to go beyond the renormalized RPA and to include vertex renormalizations omitted 
in the approximation of Suhl. The simple ring and ladder approximations renormalize only the self-energy 
and use the unrenormalized interaction strength U in scattering processes. It is clear that at intermediate 
and strong coupling the actual interaction, the quasiparticles effectively feel in their mutual scattering 
events, must be renormalized by the presence of other quasiparticle pairs. Actually we have to replace 
the bare interaction /7 by a renormalized two-particle function /C|| from ( p^ projected to the interaction 
channel. We will do it systematically in the next section. 

However, there is a more basic reason to introduce vertex renormalization whenever we use mass 
renormalization. There is a relation between a two-particle function and a derivative of the one-particle 
propagator having the structure of a Ward identity. Baym nad Kadanoff formulated necessary conditions 
for a many-body theory to fulfill mass conservation. A generalized Ward identity connecting an irreducible 
vertex function with a variation of the self-energy with respect to an external potential was thereby derived 

. Interacting electron system contains not only mass but also charge. The electron-electron interaction 
represents electrostatic energy with U ~ e? ja* . In closed systems the electrostatic potential is generated 
entirely by the actual charge distribution. The charge is carried exclusively by particles involved. Hence 
a redistribution of mass density must be accompanied by a corresponding charge redistribution in order 
not to generate spurious sources of the electrostatic potential. This "charge conservation" , i.e. the entire 
electrostatic potential is generated from the charge distribution only, can be expressed for the Hubbard 
interaction as a "Ward identity" 



a»(t/,Mi.) 
dU 



E 



E 



T 



[Kii - Xii] +"iT»^ii 



(27a) 



with grand potential fl defined in and Kii, Xii as the static, local compressibility and susceptibility, re- 
spectively. We can generalize this static identity using small space and time inhomogeneous perturbations 

U ^ U + SUij{T,T') and IJa ^ ^J■a + Smcr{T) tO 



(5$[C7,G] 



5C/ij(r,r') 



su=o 



su=o 

Sfj. = 



(27b) 
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Both these identities are fulfiUed for an exact solution. The left-hand side of ( 27b ) is the fundamental 
two-p article function for which the explicit linked-cluster expansion is used. The right-hand side of 
(27b) makes connection to the dynamical susceptibility and compressibility obtained as a variation of 
the one-particle Green function. If (27b) is violated we cannot refer to inequality (0). Identity (27b) is 



important for the linked-cluster expansion where it is used as a fundamental tool for the construction 
of the generating functional, i.e. for explicit integration in the linked-cluster theorem (^) js^. However, 
even for classical spin systems, where the situation is much easier, it was impossible to integrate ( |27b| ) in 
full and to resolve the generating functional $ completely Q . 

Most of the existing self-consistent approximations renormalize only mass, i.e. use renormalized one- 
electron propagators neglecting vertex renormalization. The Ward identity (^^ is always violated in 
self-consistent theories that do not renormalize interaction strength. E.g. for the Hartree approximation 
we obtain after a Fourier transform 



^^Hartree 



(q, 



= [/X|(q, w„OXx(q,ii/™) [l - C/^X|(q,w„0^i(q,iJ^m)] 



(28) 



with the Hartree one-electron propagators. Since = Q, identity ( |27b|) is violated by terms 

proportional to U, the small expansion parameter of the Hartree approximation. Violation of "charge 
conservation" in the Hartree approximation has no qualitative impact onto the physics of the solution 
unless there is a phase transition making the right-hand side of ( p8| ) diverge. We hence cannot rely upon 
the Hartree approximation at critical points. The Hartree phase diagram must be confirmed by a more 
advance d app roximation complying better with charge conservation. At least in such a manner that both 
sides of ( ^7b| ) lead to qualitatively the same phase diagram. 

To improve the Hartree approximation towards charge conservation we use the right-hand side of 
( |2^ ) for the determination of the generating functional $[[/, G] from the linked-cluster theorem. We 
end up with a generating functional of the ring diagrams investigated in Sec. [V A . Since neither this 
aj 



roximation contains verte x re normalizations, we cannot expect Ward identity (27) to be fulfilled 
. The left-hand side of (27b) for the ring diagrams is identical with the right-hand side of (|2^ ) 
where th e Ha rtree propagators are replaced by those from the ring diagrams. The new right-hand side of 
identity ( ^7b| ) must be determined from a set of integral equations and becomes a complicated function of 
irreducible vertex functions S'Ea-/SGa-' ■ The integral equation without integration variables can formally 
be written as 



1 — G|G| 



ST, 



u 



I 

(5G| 



GiGi 



SG^ 



- G|Gt 



U 



sg; 



GiGi 



U 



6Gt 



6G 



Ring 



(29) 



The vertex functions ST^/dG^' contain first ladder terms with interaction lines renormalized due to 
electron-hole polarization bubbles, cf. Fig. 3. We see that to comply with charge conservation at least 
all ladder and ring diagrams must be taken into account. We can continue iterations beyond the ring 
diagrams by correcting the two-particle function C(q, Wm) using (p9|). From the new function C(q, it'm) 
we derive the new self-energy and the vertex functions ST^ / 5Ga' . They again generate new corrections 
to C(q, and the iterations can continue till the convergence is reached. Already the initial two steps 
of this iteration scheme indicate that first the parquet diagrams can be expected to satisfy identity (27). 

The parquet algebra is very complicated and does not allow for a general self-consistent (non- 
perturbative) solution. If all simpler approximations break identity ( p7| ) we may ask how much the 
violation of charge conservation matters. It is physically inacceptable if the functions from the right and 
left-hand sides of equation ( p7| ) generate qualitatively different phase diagrams and lea d to incompatible 
spectral properties. We can trust an approximation violating (|2^) only if both sides of ( ^7b| ) are qualita- 
tively equivalent or if the next step in the iteration scheme towards charge conservation corroborates the 
results of the chosen approximation. E.g. the phase diagram of the Hartree approximation can be relied 
upon if it is confirmed by the ring diagrams. Phase transitions calculated with the ring diagrams can be 
seen as justified if the vertex functions 5Y,„ /SGa' do not contain non-integrable singularities that would 
generate new singularities in integral equation (^9|). 



11 



VII. DIPOLE APPROXIMATION 



It follows from preceding analysis that the full parquet approximation would be an ideal candidate 
for a theory interpolating between weak and strong coupling. It renormalizes the interaction strength 
to comply with charge conservation, at least qualitatively. It can be expected to possess appropriate 
asymptotic behavior in both the extreme limits. We, however, also learned that only nonperturbative 
solutions, treating singularities of two-particle functions at intermediate coupling analytically, can lead 
to numerically stable solutions for arbitrary interaction strength. Parquet diagrams were studied in 
connection with the local-moment problem as improvements of Suhl's approximation. Since no exact 
solutions to the parquet equations exist, approximations had to be used. The local ansatz chosen in 
failed to reproduce the Kondo scale at strong coupling. Neither non-self-consistent analysis of divergent 
diagrams at strong-coupling helped in the selection of the relevant classes of diagrams. To achieve 
a feasible extension of the simple approximations from Sec. ^ appropriate for the description of the 
metal-insulator transition we must reduce the full parquet algebra. Since the scatterings in the electron- 
hole and interaction channels contribute the same divergence at the critical point we must involve the 
two channels self-consistently. To do that we propose a new approximation scheme backed by physical 
reasoning selecting relevant scattering processes leading to the creation of neutral "dipole" bound pairs. 

In the critical region of the metal-insulator transition we expect that electrons and holes with opposite 
spins, i.e. singlets, form pairs of almost bound states. It means that if an electron and a hole with opposite 
spins get close to each other they tend to stay together for macroscopically long times. Probability of 
mutual scatterings of electrons and holes with opposite spins is hence much higher than other scattering 
events. The intra-pair scatterings of the electron and hole from the binding pairs are pronounced and are 
mediated by a renormalized interaction due to the presence of other almost bound electron- hole pairs. 
The pairs have the total charge zero and carry only a dipole moment. The inter-pair interaction is of 
dipole character and hence weak. There is no significant renormalization of the electrostatic potential 
in scarce interaction processes between pairs or pairs and particles. The dipole approximation hence 
systematically neglects all multiple scatterings where more than two particles take part. 

This intuitive physical picture can mathematically be formulated in terms of specific approximations 



in the general parquet algebra (12). First of all the three-particle and higher-order irreducible diagrams 
are neglected, i.e. the completely irreducible function I^i = U. Next, the triplet scatterings are neglected 
in favor of singlet scatterings. We can put in the asymptotic limit of the metal-insulator transition 

A^^{k,k';q)^0 (30a) 

where a — U, ee, eh. This approximation enables us to solve the parquet equations (^2|) in the vertical, in- 
teraction channel. The two-particle function from the vertical channel, i.e. /C|| (fc, k';q) = U {k, k'; q) 
becomes a functional of the two-particle scattering function from the horizontal channel defined as 

A^^ik,k':q)=IC'^^{k,k';q) = U + ICf^{k, k' ; q) + ICf^ik, k' : q) . (30b) 

In the critical region of the metal-insulator transition the two-particle scattering functions /C|| (fc, /c'; g) 

and IC^^{k,k';q) show divergence for q = qo and u ^ 0. This divergence is smeared out whenever we 
build convolution of these functions connected only by corresponding electron propagators. Only if we 
interject bare interaction U between any convolution of reducible functions from both horizontal and 
vertical channels the divergence survives in leading-order also in convolutions. This restriction reflects 
the physical fact that multiple scatterings of bound pairs with other pairs or fermions are negligible at the 
critical point where bound electron-hole pairs condense. Taking this into account we obtain a solution 
for /C|j^(fc, k'; q) in closed form and independent of the fermionic four-momenta fc, k' 

mk,k'-,q) = T,M,^.^)- ,_^l^^^^^^^^y (30c) 

The pair function A(q, w^) is a sum of irreducible diagrams that cannot be split into two separate 
components of the same character by cutting the interaction line. Typical diagrams included and omitted 
in the renormalized interaction V-\^i are plotted in Fig. 4. The irreducible function A(q, ivm) is connected 
with the two-particle scattering function IC^^{k' ,k";q) analogously to ( [l0| ) 
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k'k" 



To complete the approximatio n we must spec ify the two-particle function from the horizontal channel. 



It is determined from equations (12a) and (121) with appropriate irreducible functions. Since multiple 
electron-electron scatterings do not contribute divergent terms in the critical region of the metal-insulator 
transition we can neglect these scatterings completely without affecting leading-order asymptotics. We 
can write 

Aff(fc,fc';g)-ru(g) , Af^(fc,fc';g) = 0. (31) 



The parquet algebra is now complete, i.e. ( pl|) and ( |l2b| ) fully characterize /C|J'(fc', k"; q). To construct 
the generating functional we rewrite relation (O) between the two-particle scattering function /C and the 
correlation function C in our simplified parquet algebra to 

C(q,ii/„0 = — - — : . (32) 

1 - C/A(q,«t/„i) 

Using these equations and linked-cluster theorem (|^) we obtain the functional of the grand potential 

cjn,k 

-hG^(k,itj„)E<,(k,iw„)} + — r^2^e " / dX- .r^i...., ■ — r- (33) 

Eqs. (|3^)-(^3|) fully determine the thermodynamic as well as the spectral properties of the dipole approx- 
imation for arbitrary interaction. They represent an asymptotic solution of the parquet diagrams in the 
critical region of the metal-insulator transition, U Z' Uc, where electron-hole bound pairs condense. 

It is still a tremendous task to solve the above equations in the whole range of the coupling constant. 
In the weak-coupling limit it reduces to a sum of RPA and ring diagrams, i.e. to a FLEX approximation 
||3^ . In the strong-coupling limit the function r||((7) develops a pole at the Fermi energy and we can use 
the low-energy expansion as demonstrated for simple approximations in Sec. ^ In the critical region of 
the metal-insulator transition U —f Uc 'Vfe can replace in leading order the function A(UX) in the grand 
potential by A{U). We can integrate over the interaction strength to obtain the generating functional in 
closed form 

$[C/,G] = ^^e--°^ln[l-C/A[G](C/|q,z^™)] . (34) 

qm 



Note that ( |34| ) differs from the generating functional (15a) by a factor of two due to the fact that in the 
critical region we sum upequivalent contributions from the ring as well as from the ladder diagrams. The 
weak coupling limit of (Sj) does not hence reproduce leading-order term in U with the correct weight. 

The final step in our approximation in the critical region of the metal-insulator transition is the low- 
energy ansatz enabling to treat the singularity in C(q, C) analytically. At the mean-field level, where no 
momentum integration appears, we use 

A(C) = Ao + zsign(ImC)A'C . (35) 

We introduced two real and positive "mean-field" parameters Aq and A' that are to be determined from 
equations (p0|)-(^3|) with the above ansatz in the function r||((^). Note that there are no closed equations 
for Aq and A' derivable from a generating functional. It means that all physical quantities or equations 
for them must be determined from the full grand potential (|3|), (|3^). The low- frequency ansatz ( |35| ) 
must be applied only at the end stage of the evaluation or solution. It is important that contrary to 
the simple approximations the dynamical vertex renormalization A' is now a nontrivial function of the 
interaction strength. This is a minimal condition a theory must fulfill if metal-insulator transition should 
be described IEgII. 
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We are now in a position to perform quantitative analysis of the dipole approximation in the critical 
region of the metal-insulator transition. All the singularities in this approximation are generated by the 
single function r||(C). With the low- frequency ansatz (35) we get the singularities under control. It 
enables us to follow the solution from the weak-coupling regime up to the critical interaction where the 
metal-insulator takes place and to learn the critical behavior at this transition point in detail. However, a 
number of technical steps must be done before we gain the explicit strong-coupling asymptotics within the 
dipole approximation. Quantitative analysis and the results of the dipole approximation for the single- 
impurity Anderson model and the Hubbard mean-field model will hence be presented in the forthcoming 
paper. 



VIII. CONCLUSIONS 



We analyzed in this paper correlated electron systems described by single-impurity Anderson and 
Hubbard models at intermediate and strong couplings. We concentrated on the half-filled case where 
the Kodo behavior in the former and metal-insulator transition in the latter model are expected. We 
showed that in both cases a vanishing scale makes numerical solutions unstable and prevents numerical 
extrapolation of weak-coupling approximations to the strong-coupling limit. Hence neither Kondo strong- 
coupling asymptotics nor the metal-insulator transition are directly accessible from the weak-coupling 
(metallic) side. We demonstrated that the vanishing, Kondo, scale has its origin in a pole of an appropriate 
two-particle Green function. The relevant two-particle function at intermediate coupling is shown to be 
a dynamical generalization of the two-particle part of the underlying Hamiltonian 1^. The pole in this 
function signals that singlet bound electron-hole pairs condense and turn the half-filled metallic solution 
insulator. 

The vanishing Kondo scale at strong coupling makes quantitative description of this limit difficult. 
The only way how to get the instabilities connected with this small scale under control is to treat the 
Kondo scale separately from other scales determined by the Hamiltonian or the input parameters. To 
this end we proposed a perturbative scheme enabling direct access to two-particle functions. We used 
a quantum version of a linked-cluster expansion where bonds stand for interaction lines and vertices for 
loops of fermionic particles (electrons, holes). With the aid of the linked-cluster theorem we succeeded in 
constructing a generating functional for the parquet diagrams where all three and higher order irreducible 
functions are neglected. This approximation is a fully self-consistent sum of contributions from vertical 
(ring) and horizontal (ladder) channels of the full two-particle scattering matrix. 

The approximation containing the unrestricted parquet algebra does not allow for a global nonper- 
turbative solution needed in the strong-coupling limit. We analyzed simplified versions of the parquet 
approximation summing separately either the ring diagrams or the ladders of electron-hole and electron- 
electron multiple scatterings. These approximations, known as ring approximation, renormalized RPA 
and TMA, are used to demonstrate the two-particle instability at intermediate coupling. It was shown 
that the ring and electron-hole ladders lead to the same divergence in the effective quasiparticle mass 
in the Fermi-liquid regime. We used a low-frequency ansatz to factorize leading-order singularity of the 
two-particle scattering function. It enabled us to treat the vanishing Kondo scale analytically and sepa- 
rately from the other scales of the interacting system as demanded from numerical stability. The simple 
ring and ladder approximations were evaluated with the low-frequency ansatz in the strong-coupling limit 
and the same physically incorrect result for the rings and electron-hole ladders was obtained. The Kondo 
scale in both cases leads to the dependence Tk oc exp{— 7r^[/^i/^} either in the single-impurity Anderson 
model or in the Hubbard model in d = oo. 

The failure of the simple ring and ladder approximations to the two-particle scattering function to 
reproduce the Kondo scale in the single-impurity Anderson model or metal-insulator transition in the 
Hubbard model lies in the lack of vertex corrections or dynamical renormalization of the interaction 
strength. We used a very general argument, interpreted as Ward identity for charge conservation, con- 
necting variational derivatives of the generating functional with respect to the spin-dependent chemical 
potentials and the interaction strength. This identity, fulfilled in the exact solution, demands that 
whenever we renormalize electron mass (one-electron self-energy) we must adequately renormalizc charge 
(interaction strength). We showed that first the full parquet diagrams self-consistently summing ring and 
ladder diagrams can be expected to fulfill the Ward identity (P7|). 

Any reliable approximation suitable for the description of the transition from weak to strong cou- 
pling in interacting electron systems must contain vertex corrections qualitatively equivalent to those 
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embodied in the parquet algebra. We analyzed the parquet diagrams in the critical region of the metal- 
insulator transition from the metallic side and developed a dipole approximation. This approximation 
systematically neglects nonsingular contributions and keeps only leading-order divergent terms of the full 
two-particle scattering matrix from the parquet approximation. Together with the low-frequency ansatz 
a feasible approximation was achieved enabling to study the strong-coupling limit of correlated electron 
systems at a qualitatively new level. Moreover, the dipole approximation obeys Fermi-liquid asymptotics 
at weak coupling and it can be viewed as a theory of mean-field character interpolating between the weak 
and strong-coupling limits. Since no important, divergent, diagrams from the parquet diagrams were 
omitted in the dipole approximation, we may expect that it will correctly reproduce the Kondo strong- 
coupling behavior of the single- impurity Anderson model and eventually the metal insulator transition in 
the Hubbard model. Quantitative analysis of the sample models of interacting electrons with the dipole 
approximation is left to the forthcoming paper. 
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Figure Captions 



Fig. 1 Real part of the self-energy of the impurity Anderson model at half filli ng an d zero temperature 



calculated for the Lorentzian DOS within the ring diagrams from Sec. iVA. The first (sharp) 



extrema are clearly seen to approach the Fermi energy, {uj — 0), with increasing interaction. 

Fig. 2 First few diagrams contributing to the electron- hole, (a), electron-electron, (b), and interaction, 
(c), channels of the two-particle scattering function. 

Fig. 3 New vertex functions generated from functional derivatives of the self-energy from the ring di- 
agrams. The double wavy line is the interaction between spin up and spin down particles, and 
the double dashed line the same between the particles with the same spin renormali zed by ring 



diagrams. These vertex functions are not contained in the left-hand side of identity ( 27b ) of the 
ring diagrams. 

Fig. 4 Typical diagrams contributing to r^, (a), (b), where no more than two particles take part in 
multiple scatterings from the irreducible function A. Diagram (c) does not contribute to F-f^, since 
more than two fermion loops are multiply connected with the interaction line. 
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